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Abstract 

Population growth (or decay) in a country can be due to various f 
socio-economic constraints, as demonstrated in this paper. For exam¬ 
ple, sexual intercourse is banned in various religions, during Nativity 
and Lent fasting periods. Data consisting of registered daily birth 
records for very long (35,429 points) time series and many (24,947,061) 
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babies in Romania between 1905 and 2001 (97 years) is analyzed. The 
data was obtained from the 1992 and 2002 censuses, thus on persons 
alive at that time. We grouped the population into two categories 
(Eastern Orthodox and Non-Orthodox) in order to distinguish reli¬ 
gious constraints and performed extensive data analysis in a compar¬ 
ative manner for both groups. From such a long time series data 
analysis, it seems that the Lent fast has a more drastic effect than the 
Nativity fast over baby conception within the Eastern Orthodox popu¬ 
lation, thereby differently increasing the population ratio. Thereafter, 
we developed and tested econometric models where the dependent 
variable is the baby conception deduced day, while the independent 
variables are: (i) religious affiliation; (ii) Nativity and Lent fast time 
intervals; (hi) rurality; (iv) day length; (v) weekend, and (vi) a trend 
background. Our findings are concordant with other papers, proving 
differences between religious groups on conception, - although reach¬ 
ing different conclusions regarding the influence of weather on fertil¬ 
ity. The approach seems a useful hint for developing econometric-like 
models in other sociophysics prone cases. 


1 Introduction 

Babies are born and therefore conceived non-uniformly over the year. The 
world-wide trend is affected by seasonality (Quetelet, 1826; Lam and Miron, 
1991; Cancho-Candela et ah, 2007). The determinant factors producing these 
seasonal effects can be grouped (Friger et ah, 2009) into so called natural 
characteristics (latitude, weather conditions, day-length (DL)) and artihcial 
ones (demographic, economic, socio-cultural, including religion). 

The seasonal effects of natural factors have been extensively studied. The 
target populations are spread all around the world (Lam and Miron, 1994; 
Martinez-Bakker et ah, 2014; Friger et ah, 2009; Dorelien, 2013; Hubert, 
2014), but very rarely were they focused on Romania (except for Huber 
and Fieder, 2011). The findings on climate effects vary across time and 
geographical location. They have been extensively analyzed, the importance 
attributed to them ranging from high (Villerme, 1831) to as unimportant in 
e.g. some French population (Regnier-Loilier and Rohrbasser, 2011) or seen 
as acting jointly with economic development (Huber and Fieder, 2011; Seiver, 
1985) due to the reasoning that poor populations are negatively affected by 
low temperatures during winter or, on the contrary, by hot summers without 
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air conditioning. Other researchers (Macfarlane, 1974; Lam and Miron, 1994) 
stated that weather affects the seasonality of births as long as there are major 
temperature differences over the year. Other natural factors like day-length 
(Bronson, 2004; Rojansky et ah, 1992) or moon phases (Criss and Marcum, 
1981; Jongbloet, 1983) were also tested for their possible connection to birth 
seasonality. 

From the artihcial factors which were recorded as having some influence 
over the seasonality of births, let us mention the demographic ones: age, 
marital status, education, intergenerational effects (Bobak and Gjonca, 2001; 
Huber and Fieder, 2011) and economic ones: income, rurality (Kestenbaum, 
1987; Vitzthum et ah, 2009), as well as socio-cultural characteristics: religion, 
ethnicity (Friger et ah, 2009; Hubert, 2014). It turns out occurs that research 
on the influence of the artihcial factors is rare, perhaps due to some lack of 
precise data. In this context, our paper aims to answer the following question 
concerning the inhuence of religiosity on the seasonality of dates of conception 
and of birth. More precisely, we raise questions as: did the interdiction of 
sexual intercourse during fast (Lent and Nativity) periods have any effect on 
the date of conception in twentieth century Romania? Has that socio-cultual 
constraint an effect on population growth? 

2 Data sample 

Romania is a country with a large majority (more than 86% nowadays) of 
Eastern Orthodox population as can be observed in Fig. 1, which is based on 
censuses performed across Romanian territories in the last 150 years. Other 
religious denominations have been present in different proportions during the 
last 150 years in regions that are now part of Romania. For example, it is 
apparent that Transylvania is the most heterogeneous Romanian territory 
with regard to religious diversity (Herteliu, 2013). However, this effect is 
apt to be minor, and our current analysis is conducted at the national level, 
without distinguishing regions. 

It is known that the operation of complex mechanisms acting on or re¬ 
lated to religion is not an easy task to investigate (Herteliu, 2007; Herteliu 
and Isaic-Maniu, 2009; Ileanu et ah, 2012). In this context, to get some sta¬ 
tistical balance given the vast Eastern Orthodox majority, the non- Orthodox 
segments are grouped together, thus transforming the religious affiliation into 
a simple yes or no variable. Even if some other religions grouped into the 
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[ BOrthodoi OCtholic aRocnanC<thoi»c DWetorm»d oCie<fcCatfiolic lUnttirian oJtirith oLuihwan lOthw | 



Figure 1: Religious affiliation in Romania Source of data: Censuses per¬ 
formed across Romanian territories in the last 150 years. Please note that 
in communist era four censuses were performed (1948, 1956, 1966 and 1977). 
Due to political bias, none of them included the questionnaire items regarding 
religious affiliation. Note: The data on Romanian territories are organized 
or grouped differently from today. The codihcation is the following: T Tran¬ 
sylvania; P Romanian Principalities (Moldavia and Wallachia); GR Greater 
Romania 
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category non- Orthodox have periods of fasting during the year, due to their 
particularities - such as the use of a different calendar, different dates for 
Easter, different rules for fast periods, or other features specihc to various 
religious sects or denominations, their exclusion from the well dehned Eastern 
Orthodox group was preferred for the sake of simplicity. 

It is rather well-known that in the Eastern Orthodox tradition, in ad¬ 
dition to abstaining from particular types of food and drinks during fast¬ 
ing periods, individuals should also avoid sexual intercourse. Support for 
these extensive fasting interdictions specihc to the Eastern Orthodox tradi¬ 
tion {http : / / orthodoxinfo.com/praxis/ fasting sex.aspx retrieved in March 
2015) are based on interpretations of certain quotes from the Bible (such as 
Romans 8:12-14 or Galatians 5:16- 17). One of our goals is to hnd whether 
these constraints have been implemented. 

2.1 Data sources 

The analysis is based on a very long (35,429 points = 365 (days/year) ? 97 
years -|- 24 leap years) daily time series representing all births (24, 947,061) 
of persons alive at the 1992 or 2002 censuses (similar to Kestenbaum, 1987). 
Other long data series have been studied before (Cancho-Candela et ah, 2007; 
Ausloos et ah, 2015; Rotundo et ah, 2015) using different methodologies. In 
the present case, data points were recorded for 01/01/1905 to 31/12/1991 
(31,776 points) from the 1992 census, while the other 3,653 points (corre¬ 
sponding to the period 01/01/1992 to 31/12/2001) are from the 2002 census. 
The data source is the Romanian National Institute of Statistics (NIS) via 
a query tool available within NIS’s intranet [http : //happy : '&1/PHC) re¬ 
garding 2002 and 1992 censuses. For convenience, in the current research, 
the terms born per day” will be used instead of ’’number of births on a 
given day for persons still alive at the 1992 and respectively at 2002 Cen¬ 
sus”. The data for DL and moon phase (proportion of moon disc illumina¬ 
tion) was obtained from http://www.timeanddate.com, For daily tempera¬ 
tures the source for data is Romanian National Meteorological Agency (via 
http : //www.ecad.eu/dailydata/customquery.php (retrieved on December 
2013)). The dates for Easter in the Eastern Orthodox religion during the last 
110 years is obtained from http : //www.smart.net (retrieved on December 
2013). 
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2.2 Variables 


(a) Variables from censuses : Data was based on persons’ self-declared birth¬ 
day. In addition, we use data about: (i) religion: distinguishing Eastern Or¬ 
thodox (EOx) and Non-Orthodox (NOx) affiliation. We must assume that 
the religious affiliation of the parents is well dehned, based on the informa¬ 
tion given by their offspring at census time, and is the same for both par¬ 
ents. Even so, the available statistical information showed that the intensity 
of changing religious affiliation across the Romanian population is very low 
(Ileanu et ah, 2012); (ii) rurality: the percentage of people located in rural 
or non rural areas. This introduces a double filter: Eastern Orthodox (rural¬ 
ity EOx: rEOx) and Non-Orthodox (rurality NOx: rNOx), the alternative 
being Eastern Orthodox in Urban communities (EOxU) and Non-Orthodox 
in Urban communities (NOxU). 

After subtracting 280 days from the registered birthday, an estimate of the 
conception day was obtained (Seiver, 1985; Lam et ah, 1996). However, date- 
of-birth evidence can only point to to sexual intercourse leading to successful 
pregnancies. It is known that Lent is 48 days long while the Nativity fast 
lasts 40 days. The uncertain conception date, since not every child is born 
exactly 280 days after conception (Pana et ah, 2015), overlaps the 40 or 
48 day fast. We assume that the birth rate is correlated with the rate of 
intercourse. Practically, family planning and contraception methods have 
not been very common in Romania (it was certainly not encouraged by the 
communist regime between 1966 and 1989) (Chelbegean, 2010). We assume 
that the birth rate is correlated with the rate of intercourse. Recall that 
this a posteriori estimate of sexual intercourse resulting in a birth a is made 
primarily to correlate to natural variables, i.e., weather conditions and day 
light (or night) duration at the time of the sexual activity. 

For each year, the fraction of births yi on given day, i, were computed 
separately for the orthodox , and non-orthodox populations y^^^, using: 


yr^ = 


n. 


EOx 




EOx 

i 


( 1 ) 


and 


Vi 


NOx 


^NOx 




NOx 

i 


( 2 ) 


where and is the number of births on the specihc day i for EOx 

and NOx respectively, while is the total number of 
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births during that year for the respective populations. 

The dependent variables in the regression models are calculated by divid¬ 
ing the daily fraction of births (1) and (2) by the daily expected (uniform) 
distribution (UD) across the year which is = 1/365 for non-leap years 
and = 1/366 for leap years. Therefore, we dehne the following two 
dependent variables: 




,EOx 


y 


UD 


(3) 


and 






NOx 


,UD 


(4) 


respectively and using the appropriate y^^. 

(b) Nativity and Lent fast variables: We used a countdown variable (in 
days) to Christmas or Easter for each day of Nativity or Lent fast, instead 
of a simple binary variable indicating whether or not conception occurred 
in the fasting period. We preferred this because it is possible that religious 
constraints were observed more scrupulously the closer one was to either 
Easter or Christmas. This countdown variable was computed having the 
conception date as the reference point. Subsequently, Days Before Christmas 
(DBC) was computed from: 


DBC = 0 or t (5) 

if the conception occurs within or outside Nativity fast time interval; t rep¬ 
resents the number of days to Christmas. 

Similarly, Days Before Easter (DBE) was dehned but where t represents 
the number of days to Easter: 


DBE = 0 or t (6) 

For completeness, we note that Romania adopted the Gregorian calendar 
in 1919. Therefore, the 31st of March 1919 was followed by the 14th of 
April 1919. Moreover, the Romanian Eastern Orthodox Church is using an 
updated form of the Julian calendar. Thus, Christmas is celebrated in Ro¬ 
mania on the 25th of December. However, Easter for the Romanian Eastern 
Orthodox Church is usually not synchronized to that of the Roman Catholic 
Church. Even so, it is emphasized that these specihc calendar change details 
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do not interfere with the results of our research, since the data regarding reg¬ 
istered birthdays in the 1992 and 2002 censuses were based on the Gregorian 
calendar. 

(c) Time interval related variables : Some binary variables were inserted 
in order to mark some particular time windows during which the concep¬ 
tion date might have occurred and influenced by global events or days of the 
week. These are: (i) First World War (IstWW) taking the value 1 for con¬ 
ception days between 28 July 1914 and 11 November 1918 and 0 outside this 
timespan; (ii) Second World War (2ndWW) having a 1 value for conception 
days between 1 September 1939 and 2 September 1945 and 0 outside this 
timespan; (iii) weekend’s days having a value = 1 if the conception day is 
on a Saturday or Sunday and 0 for Monday, Tuesday, Wednesday, Thursday 
and Friday and (iv) squared trend for t=l to 35,429). 

(d) Natural variables : The effect of some natural variables is also ana¬ 
lyzed through: (i) day length (DL) measured as the number of minutes of 
daylight during each 24 hours of a day. This information refers to Bucharest 
as the capital city of Romania, but is taken as to be the same for wher¬ 
ever place in Romania; (ii) temperature in Celsius degrees (mainly mea¬ 
sured for Bucharest) (Due to the usual intrinsie error bars on meteorologieal 
measurements (and variation aeross the country and throughout the day), 
using Bucharest temperatures is just a rough approximation to local temper¬ 
ature.), and (iii) proportion of the moon disc being illuminated (having also 
Bucharest as the reference point), within the same approximation. 

2.3 Statistical analyses 

The analysis began in qualitatively to determine whether births (or concep¬ 
tion) vary throughout the year. Graphs were drawn with various softwares 
(Kaleidagraph, EViews, Excel, Power Point). Subsequently, the daily av¬ 
erage conception within and outside fasting periods was computed. Prior 
to the econometric modeling, we checked a possible periodicity via a Fast 
Fourier Transformation performed with SPSS for 2^^ cases (the last 32,768 
days). 

The hnal step in the data analysis consists in the design of several econo¬ 
metric models. Breaking points in the trend of the dependent variables were 
identihed with Cusum SQ Test performed with EViews. Regressions were 
performed with IBM-SPSS (Ordinary Least Square - OLS estimations), - the 
same model being run for each phase/cycle (identihed after a Cusum SQ 



Indicator 

EOx 

NOx 

N Valid 

366 

366 

Missing 

0 

0 

Mean 

59219.15 

9082.180 

Median 

59418.00 

9050.0000 

Std.Deviation 

8906.85 

1049.23 

Skewness 

0.93 

2.21 

Kurtosis 

13.52 

37.59 

Minimum 

11291 

1904 

Maximum 

131311 

20016 

Sum 

21629796 

3317625 

Percentiles 20 

53554.80 

8452.8 

40 

57996.20 

8851.8 

60 

60449.60 

9209.2 

80 

65141.40 

9.618.6 


Table 1: Descriptive statistics for daily aggregated data of population growth 


Test). At this stage, factors suspected of being collinear were removed (e.g. 
day length was chosen instead of temperature) from the regression models. 
The same procedure was performed for factors with low statistical signih- 
cance (e.g. moon phases). 

3 Results 

3.1 Preliminary data inspection 

As shown by Huber and Fieder (2011) for women born between 1920 and 
1955 and older than 45 years (completed fertility), Romanian births register 
a non-uniform distribntion across the year with a maximnm in Jnne and 
minima in December and January, as obtained from monthly data. 

Onr raw data contains the number of births on each day of the year (EOx 
and NOx) for a long period (97 years). In order to present this extensive data 
in a concise way, a pivot table was designed separately containing monthly 
and daily aggregated valnes for EOx and NOx. For both populations, the 
mean and median are very close to each other (less than 0.4%) while the stan- 
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dard deviations for the means are quite low (15% of the mean for EOx and 
11.6% for NOx). The statistics of both populations are slightly asymmetric 
and leptokurtic. Other descriptive statistical characteristics are available in 
Table Next, the aggregated date were sorted into quintiles (Dedu et ah, 
2014; Giuclea and Popescu, 2009) obtained by quintiles. The following inter¬ 
vals were obtained for EOx (Fig. |^; (ql) days having under 53,554.8 births 
labeled in red; (q2) days having between 53,554.8 and 57,996.2 births labeled 
in orange; (q3) days having between 57,996.2 and 60,449.6 births labeled in 
yellow; (q4) days having between 60,449.6 and 65,141.4 births labeled in light 
green and (q5) days having over 65,141.4 births labeled in dark green. Sim¬ 
ilarly for NOx (Fig. |^, the quintile intervals are (ql) days having under 
8,452.8 births (red); (q2) days having between 8,452.8 and 8,851.8 births 
(orange); (q3) days having between 8,851.8 and 9,209.2 births (yellow); (q4) 
days having between 9,209.2 and 9,618.6 births (light green) and (q5) days 
having over 9 618.6 births (dark green). 

This sorting by quintiles of almost one century’s data can (by avoiding 
compensation through averaging over unequal size month) provide some sta¬ 
tistically signihcant differences among daily distributions across the year. It 
can easily be seen that this compensation does not presently occur is our 
case. Therefore, it may be deduced that the outlined monthly or daily dif¬ 
ferentiation does not lead to a random outcome. 

Referring to Fig 2, for the EOx population, there are important differ¬ 
ences between November, December and January and the rest of the year, 
and also between the hrst two thirds of a month compared to the last third. 
These monthly differences were tested by a test and are statistically sig¬ 
nihcant (computed was 131,285.4 being signihcant with a p value lower 
than 0.001). Again referring to Fig. 2, the lowest aggregate number of births 
for a specihc day was, naturally, on 29th of February (11,291 births) while 
the maximum was on the 1st of January (131,311 births). We will comment 
on the large number of births on 1 Jan below. 

In the case of the NOx population, (See Fig. 3) the lower birth hgures 
in November and December are consistent with lower number of births in 
those months in the EOx population, while January is not so distinct as in 
the NOx case. The distribution of the daily hgures outside the November- 
January time span seems to be more equally balanced for NOx, - no persistent 
pattern is visible, in contrast to the EOx case. Similarly in the EOx and the 
NOx population, the minimum aggregate number of births was registered on 
29th of February (1,904 births) while the higher one occurs on 1st of January 
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Model 

Int 

ere 

ept 

Dec 

DBE 

DL 

Rurality 

(ORTHD) 

Trefid*2 

Week-end 

Istww 

2MWW 

Adjuste 

dR 

Square 

Model 

validity 

(SigF) 

MOOveral 

190S-2001 


0.002 

-0332 

0.051 

0.073 

0.087 

0.028 

■0.055 

•0.004 

0.127 

0.00 

Ml 

190S-1910 


•0.0849"’ 

-0J676’” 

-0.0596'" 

03414’" 

0.1557’ 

0.0344"’ 



0.187 

030 

M2 

191M9U 


•0.0360 

•0.2534 

0.1252 

03107 

0.0586 

0.0838 

•0.0715 


0.229 

0.00 

M3 

1919-1947 


0.0194" 

-03034’" 

0.0105 

0.8986"’ 

0.0578"’ 

03427"' 


-0.048 " 

0.323 

030 

M4 

1949-1966 


0.0271 

•0.3457 

0.0276" 

0.1139 

0.0510 

■0.0574 



0.146 

030 

MS 

1967-1989 


0.0487"’ 

-03897" 

0.1684"’ 

0.1020"’ 

■0.0313’" 

0.0445"’ 



0.127 

030 

M6 

1990-2001 


0.0577’" 

•03316’” 

03172"’ 

0.0653"’ 

■03841’" 

0.0207"’ 



0335 

030 


Table 2: The logarithm of the ratio of conceived orthodox babies relative 
to assumed uniform daily value for 1 year. The regression coefficients shown 
are the standardized ones (prior to estimation all independent variables are 
standardized). Factors (independent variables) are: Days Before Christmas 
(DBC), Days Before Easter (DBE), Day Length in minutes (DL), Rurality 
for ORTHD, squared trend {t^ for f = 1 to 35,429), Binary variables: (i) 
FirstWorldWar (IstWW)taking the value 1 for conception days between 28 
July 1914 and 11 November 1918 and 0 outside this timespan; (ii) Second 
WorldWar (2ndWW) having a 1 value for conception days between 1 Septem¬ 
ber 1939 and 2 September 1945 and 0 outside this timespan; (iii) weekend 
days having a value 1 if the conception day is on a Saturday or Sunday and 0 
for Monday, Tuesday,Wednesday, Thursday and Friday. Bold values denote 
the models goodness of fit, adjusted by number of independent variables. 
* Statistically signihcant at level 10%; ** Statistically signihcant at level 5%; 
***Statistically signihcant at level 1%. 
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Model 

Inter- 

CM 

D6C 

D6E 

DL 

(NORTHD) 

Trend-'! 

Weekend 

IstWW 

2ndWW 

Adjusted 

R Square 

Model 

validity 

(SigFI 

MO Overall 

1905-2001 


0.007*** 

-0.157*** 

0.0071** 

0.006*** 

0.0020*** 

0.042*** 

-0.094*** 

-0.004 

0.047 

0.00 

Ml 

1905-1910 


0,0154 

-0.0420** 

■0.0285 

0.0980"* 

0.0252 

0.1007'" 



0.023 

0,00 

M2 

1911-1918 


0,0438" 

■0.1401"' 

0.0942"' 

0.1031"* 

0.1143*" 

0.0953'" 

■0.2754'" 


0.092 

0,00 

M3 

1919-1947 


0.0831"* 

■0.1794*" 

0,0304'" 

0.1712"* 

0.0890*" 

0.1112'" 


0,0010 

0.073 

0.00 

M4 

1948-1966 


0.0658'" 

■0.1919'" 

0.1000'" 

0.0490"* 

-0.0338'*' 

0.0688'" 



0.055 

0.00 

M5 

1967-1989 


0.1052'" 

-0.1888'" 

0,1948"' 

•0.0035 

-0,0023 

-0.0517'" 



0.072 

0,00 

M6 

1990-2001 


0.0969'" 

■0.1783"' 

0.1040"' 

0.0013 

0.0240 

■0.2468"* 



0.10S 

0,00 


Table 3: The logarithm of the ratio of conceived Non-Orthodox babies 
relative to assumed uniform daily value for 1 year. The regression coeffi¬ 
cients shown are the standardized ones (prior to estimation all independent 
variables are standardized). Factors (independent variables) are: Days Be¬ 
fore Christmas (DBG), Days Before Easter (DBE), Day Length in minutes 
(DL), Rurality forNORTHD, squared trend (f^ for f = 1 to 35,429), Binary 
variables: (i) FirstWorldWar (IstWW)taking the value 1 for conception days 
between 28 July 1914 and 11 November 1918 and 0 outside this timespan; (ii) 
Second WorldWar (2ndWW) having a 1 value for conception days between 
1 September 1939 and 2 September 1945 and 0 outside this timespan; (iii) 
weekend days having a value = 1 if the conception day is on a Saturday or 
Sunday and 0 for Monday, Tuesday,Wednesday, Thursday and Friday. Bold 
values denote the models goodness of ht, adjusted by number of independent 
variables. **Statistically signihcant at level 5%; *** Statistically signihcant 
at level 1%. 
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Legend 


{Under 53 554.8 
1 60 449.6-65 141.4 


1 53 554.8-57 996.2 
lover 65 141.4 


57 996.2-60 449.6 


Figure 2: Quintile distribution of daily aggregated births for the Eastern 
Orthodox (EOx) population during 1905-2001. 


13 
































































Legend 


|Under 8452.8 
■ 9209-2-9618.6 


1 8452.8-8851.8 
lOver 9618.6 


8851.8-9209.2 


Figure 3: Quintile distribution of daily aggregated births for the Non- 

Orthodox (NOx) population during 1905-2001. 
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(20 016 births). For this group, the computed 'w^is 9,184.2 showing that 
the differences are statistically significant with a p value lower than 0.001. 

For both (EOx and NOx) populations. Figs 2- 3 show a noticeably high 
number of births (almost dark green) on the 1st day of any month and a low 
number of births) on the last day of each month. 

As noted above for both populations, there is an anomalously high num¬ 
ber of births on Jan. 1st. In fact there are twice as many births on Jan 1st 
than on Jan 2nd and almost 4 times as many births as on the 31st of Decem¬ 
ber. We call this anomaly the 1st of January effect; further work should be 
done on this subject. A similar effect was noted by MacFarlane (1974). In 
brief, possible explanations for this effect might be: (i) parents’ psychologi¬ 
cal comforting thought that their child is considered to be one year younger 
if registered on January 1st instead of December 31st (a quite comfortable 
situation for both girls and boys since there would be a delay of one year to 
the compulsory military training, - effective in Romania prior to 2007); (ii) 
municipal recorders not working during winter holidays therefore increasing 
the likelihood of an incorrect record (despite the fact that children can be 
registered on the 5th of January mentioning the correct birth date); (iii) 
on the 1st of January there is a very important holiday for the EOx (Saint 
Greater Basil), evidenced by many children born around this day named 
Basil. (According to http : //www.capital.ro/peste — 600000 — romani — 
isi — serbeaza — ziua — onomasticade — sfantul — vasile — 142714.htm/, 
retrieved on December 2014, there are almost 600 000 people in Romania 
named Basil or the female version); and (iv) mistaken registration. (Even 
if it is not statistically important, one co-author of the current paper has 
a brother registered as born on January 1st, even if the event certainly oc¬ 
curred, on 2nd January 1:00 AM). 

Since the variation of the number of births across the year is undoubtable, 
we now address the primary concern of this research: did major fast periods 
(i.e. before Christmas and Easter) affect conception or not ? First, note 
that several religions (Roman Catholic, Greek Catholic, Old Rite Eastern 
Orthodox) grouped within NOx category still have different Lent and Nativ¬ 
ity fast periods. Moreover, the fast periods (especially for Lent) are usually 
not similar to EOx. Therefore, a visual inspection of this variability (Figs. 4 
- 5) was conducted, only for EOx. Prior to this, it seems useful to compute 
for each year the daily average conceptions within such fast (or out of fast) 
periods. 

Despite the long time series which certainly must depend on many local 
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Figure 4: Daily average number of conceptions during Lent fast/out of fast 
period for Eastern Orthodox (EOx) 


or general factors, Fig. 4 shows that the daily mean number of conceived 
children by orthodox is significantly lower during Lent than during the rest 
of the year. In the whole range of almost 100 years only two exceptions to 
this rule can be seen: (i) during the First World War (more exactly, in 1916) 
and (ii) during 1989, a year which had a strong emotional impact on the 
Romanian population, when there was a revolution and a major change of 
political system. Even if births increased because of Ceausescu’s decree, the 
number of conceptions out of fast is still larger than the number of concep¬ 
tions during fasting periods. Thus, even if in the short term (1967-1968) 
the Ceausescu Decree reduced the gap between ’’out of fast” and ’’during 
the fast” on the long term, the decree had not enough force to completely 
eliminate the effect of ecclesiastical admonitions. 

Similar plots of conception within and without the Nativity fast (Fig. 5) 
show dramatically different behavior. In this case the rule of fewer concep¬ 
tions during fast periods was broken in 67 from all 97 available years. To 
complete the basic statistical analysis, an enhanced econometric analysis is 
described and performed in the Econometric models sub-section, below. 
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Figure 5: Daily average number of conceptions during Nativity fast/out of 
fast period for Eastern Orthodox (EOx) 

3.2 Identification of periodicity 

The analytical approach of our paper is an econometrical one. However, 
before applying econometric model, it is useful to test whether the main 
dependent variable of interest EOx is affected by some periodicity (figure 
6) in order to avoid such a specific trend. A spectral analysis is thereby 
performed, using SPSS (Tukey-Hamming method Brillinger, 2002). Such a 
spectral analysis may be done when the number of data points is a power of 2. 
Therefore, in order to have the maximum set of points, 2^^ (= 32,768 cases), 
the first 2 661 data points were deleted, whence shortening our dataset to 
the 14/4/1912 - 31/12/2001 interval. 

The major peak (we use a log 2 scale in order to slightly enhance the x- 
axis) is about half a year (183 days or a little bit more than 26 weeks). The 
next largest peaks occur ca. 372 days (ca. 53 weeks), four months (120 days) 
and one month (30 days); such ’’periods” are quite similar to those found in 
(Cancho-Candela et ah, 2007). 
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3.3 Econometric modelling 

As presented in the above paragraphs, the number of conceptions in a pe¬ 
riod of time was assumed to be determined by a set of factors, more or less 
known and more or less understood. In the current analysis, classical multi¬ 
ple regression models are considered, with various components and variables 
measured on a daily time scale. Several versions of the models were tested 
starting from the general linearized multifactorial model (Andrei and Bour- 
bonnais, 2008). All factors mentioned in the data and methods section above 
were initially included in the models. Based on the goodness of fit of these, 
taking into account the statistical signihcance of the regression coefficients 
(with a Student t-test), an optimum, model was obtained and it will only be 
the only model outlined below. As in other research (e.g., Lam and Miron, 
1996), the best results are obtained when considering the logarithmic variable 
as a dependent variable. 

(A) Model for the orthodox population; components and variables: 

In such a case, the dependent variable, and covariates are: Days be¬ 
fore Christmas (DBG), Days before Easter (DBE), rurality EOx (rEOx), day 
length (DL), World wars, through IstWW and 2ndWW variables, weekends 
(WE), trends (long term tendencies, a stable component, core of the time 
series) and a stochastic component (e). The variables are described in a sub¬ 
section of the Material and methods section above. A formalized version of 
the model can be easily written; after applying a logarithmic transformation, 
the model reads in obvious notations, ln{d^^^) is equal to 

= ttg T DBC -\- (y .2 DBE -\- 03 ln(^D 

-|-q ;4 rEOx + 0:5 IstWW + ae 2ndWW + aj WE + ag {trendy + e (7) 


(B) Model for the non-orthodox population; components and variables: 
with almost the same explanatory variables but instead having the depen¬ 
dent variable, d^^^, i.e. the normalized number of non-orthodox conceptions. 
Moreover, instead of covariates share of orthodox in the rural area, the share 
of non-orthodox in the rural area is introduced (rNOx), which may be more 
suggestive. Thus, in obvious notations, ln{d^^^) is equal to 

= Po + /3i DBC + P 2 DBE + /^s ln{DL) 

+/34 rNOx + 135 IstWW + /Sq 2ndWW + WE + /dg {trendf + e (8) 
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Figure 6: Fast Fourier Transformation (FFT) for the This dependent 

variable in the regression models is calculated by dividing the daily fraction 
of births by the daily expected (uniform) distribution across the year. 
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Figure 7: Breaking points identification using CUSUM SQ test 


(C) Structural break identification: 

Since the series covers a very long time interval (97 years=35,429 days), 
parameters might vary somewhat. Perhaps as a result the overall model (MO, 
presented in Table 1) appears not to be very efficient in that some variables 
such as the Christmas fast seem to be insignihcant. As in (Cancho-Candela 
et ah, 2007), a segmentation of the long time span into shorter sub-periods 
seems an interesting procedural test. 

In order to do so, the CUSUM SQ (Young, 2011) test was applied to the 
overall estimated model (MO). The CUSUM SQ graphic, presented in Figure 
7 suggests that a 5% level of signihcance multiple changes is likely during the 
analyzed timespan. It can be observed that the breaking points are ca. 1911, 
1918; 1947, 1966, 1990. Consequently, there are 6 suggested sub-periods of 
stability for the models (Mi): 1. 1905-1910 (Ml); 2. 1911- 1918 (M2); 3. 
1919-1946 (M3); 4. 1947-1966 (M4); 5. 1967-1989 (M5); 6. 1990-2001 (M6). 
Each break point is in fact well related to a social/ political/historical event: 
e.g., during the period for 1907-1910, there were peasants’ uprisings; at the 
end of the second and third sub-period, there were events which determined 
the decrease of the number of births: the First World War and the Second 
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World War respectively. Also, this observation suggests how to specify the 
trend component of the models. Therefore, the econometric models were re¬ 
worked for each of the sub-periods taking into account these remarks: (i) for 
model Ml a non-linear (let us take it squared, within usual series expansion 
approximations) trend of births during the entire period is introduced; (ii) for 
model M2 and M3, similar squared trends can be also assumed.; (hi) the same 
assumptions are introduced for models M4 and M5. According to national 
reports, during the period of 1956-1966, the nativity rate decreased dramat¬ 
ically due to many factors such as: abortion liberalization, better access to 
work and education for women; moreover, between 1985 and 1990, the social 
condition of Romanians signihcantly decreased. The Decree of Ceausescu 
outlawing abortion (1st October 1966) had effects only for a short period 
of time; (iv) dually, for the 6th model, the trend in the number of babies 
conceived and born decreased according to many factors, like abortion liber¬ 
alization (992 000 such operations occur in 1990), but also post-communist 
socio-economic crisis, migration etc. 

After applying estimations using SPSS methodology, the optimized re¬ 
sults are given in Table 2 and Table 3. 

After performing a F test - ANOVA it has been found that the model 
for every sub-period for EOx and NOx population are valid and statistically 
significant. The proportion of variance (R2) explained for the dependent 
variable is between 12.7% (MO and M5) and 32.3% (M3), - values similar to 
other studies (Friger et ah, 2009). Almost all regression parameters are sta¬ 
tistically significant (from a Student t-test) with slightly different degrees of 
significance (most of them having the p value lower than 0.01). As expected, 
the model performed on NOx population, registered levels of R2 lower than 
the previous ones. The maximum deviation occurs for M6 (10.5%) and the 
minimum for M2 (2.3%). Regression coefficients in this case fail to be slightly 
more often signihcant. 

In addition, the sensitivity of the statistical models was tested by taking 
into consideration a gestation duration shorter by one and two weeks, respec¬ 
tively. The models outcomes in these scenarios are found to be quite similar 
to the presented ones. 
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4 Discussion and Conclusions 


Romanians following Eastern Orthodox religious beliefs or rules, i.e. the ma¬ 
jority of religious persons in Romania, are taught to avoid sexual intercourse 
and to abstain from particular food and drinks during fasting periods. The 
results of our time series analysis on birth rates show different behaviors 
during the major fasts. Conceptions, during Lent Fast, are consistent with 
religious constraints. The negative value of the coefficient (DBE) for Days to 
Easter during the Fast shows that as the Holiday gets closer, the estimated 
number of conceptions signihcantly decreases. It is worth mentioning here 
that the Orthodox religion grants some exception periods to believers for rea¬ 
sons mainly linked to health. In these grace periods, there might be increased 
conception. It is important to note that the standard absolute value of the 
coefficient of the variable which measures the impact of Lent fast on con¬ 
ception (DBE) remains relatively stable in [-0.25; -0.35] for all sub-periods. 
This points to the continuity and persistence of the Eastern Orthodox faith 
in general, and particularly regarding interdictions of Lent fasting through¬ 
out the almost 100 years of this study. This continuous adhesion to faith 
is also highlighted by the values estimated in the second model (Table II) 
for the Non-Orthodox sub sample (NOx) where the estimated coefficient of 
the variable has reasonable values, though half [-0.04;-0.19] that for the EOx 
population. On the other hand, the Nativity Fast no-sex constraint does 
not seem to have a great signihcance for Orthodox Romanians nowadays: 
the coefficient of the variable (DBG) which measures this impact is positive 
after the 1st World War, showing that as the Christmas day gets closer, the 
number of conceptions increases signihcantly. If in the case of the Lent Fast 
the coefficient remains stable in time, in the case of the Nativity Fast, the 
standard value is increasing in time, showing that when approaching the Na¬ 
tivity date, the population practices a more liberal lifestyle, obeying fewer 
constraints during the fasting period. For example, in the period 1990-2001 
the value of the coefficient is approximately three times higher than the one 
estimated for the period before the Second World War. This can be likely 
explained from a rural habit: weddings are traditionally performed in ru¬ 
ral Romania in the last part of autumn, after the harvest (Trebici, 1979). 
Moreover, the fermentation process of the new wine is stopped, thus render¬ 
ing the liquid more drinkable, in particular at weddings. It has been shown 
(Trebici, 1979) that usually about 9-10 months later, the hrst baby of a new 
couple is delivered. Another point is based on the psychological idea that 
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the Christmas holiday is a more joyful event compared to Easter (in agree¬ 
ment with others, like Seive (1985)). We can also propose as a cause that, 
in rural Romania, there are fewer household duties (see Lam and Miron, 
1991; tangential Becker, 1991) occurring outside the house during the Nativ¬ 
ity fast than during spring Lent. And lastly, the day length (and of course 
night length) needs to be mentioned as there was a lack of electricity until 
the middle of the 20th century and also a lack of other indoor sources of 
entertainment, thus increasing the likelihood of sexual encounters at those 
times. 

Therefore, the proposed models highlight the impact of religiously con¬ 
trolled fasting periods on conception, besides the existence of other factors 
with persistent or temporary influence. For instance, the day length variable 
reveals the signihcant influence of several factors such as: the regularity of 
daily activities determined by the season and the temperature (these two 
factors were included in the initial models and eliminated due to the multi¬ 
colinearity effect). With one exception (Ml), our analysis reveals a photope¬ 
riodicity related cause for baby conception as found for European Countries 
like France, England, Sweden (Lam and Miron, 1994) on monthly data. Of 
course, due to multicolinearity, the photoperiodicity could have only an ap¬ 
parent effect as a proxy for factors like the temperature, outside household 
duties or the lack of available inside sources of entertainment. 

Weekend is a factor inserted into our models that tends to function as 
proxy for holidays, and more permissive sexual activity. Previous research 
(Bobak and Gjonca, 2001) suggests that during holidays, the frequency of 
conception seems to be higher indeed. Our data analysis generally agrees 
with these Endings. Therefore, for EOx population, weekends seem to have 
a positive impact on conception with the exception of M4. Furthermore, the 
same influence is reflected on NOx population also until 1966 (models Ml to 
M4). Thereafter, weekends are shown to have a negative influence. 

It also needs to be stated that despite other researchers’ Endings (Jong- 
bloet, 1983), the proportion of illuminated moon disc was found to be in¬ 
significant from a statistical view point (along a Student t-test). 

Furthermore, the values of the regression coefficients specific to the Or¬ 
thodox persons in the rural environment, illustrate that behavioral patterns 
in these geographical areas are much more prone to seasonal elements in 
child conception as compared to urban areas, where these differentiations 
are diminished. 

The trend of the series for each analyzed sub-period confirms the expec- 
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tations stated above. Thus, in each major timespan here analyzed, there are 
years of maximum increases followed by a descending evolution. The only 
exception to this rule are the years between 1948-1966, when the maximum 
points of the number of conceived babies are at the ends of the interval. On 
the one hand, this could be determined by the necessity to cover the disaster 
caused by the war; on the other hand, it shows the signihcant impact of 
Ceausescu’s decree. 

Finally, the splitting of the populations into EOx and NOx categories is 
proven relevant. The model for all sub-periods analyzed is much weaker (from 
econometrics viewpoint) for the NOx population which does not have Nativ¬ 
ity or Lent fasts, or if it does, with time delays versus the EOx population. 
Moreover, the standardized coefficient that measures the impact of Lent Fast 
on conception has in general the largest absolute value. It is followed in ab¬ 
solute value by the rurality (rEOx) coefficient, both values strengthening the 
value of the orthodox tradition and faith, factors which are underestimated 
or missing in previous human conception researches. 

In brief, the primary result of the current investigation is that the non¬ 
uniformity of the conception of babies, in Romania, for the last century, 
across the year is affected not merely by seasonality, but by a double hlter: 
season and religion. 
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